7rAE Coupling Extracted from Hyperonic Atoms 
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The latest measurements of the atomic level width in E-hyperonic Pb atom offer the most accurate 
datum in the region of low-energy E-hyperon physics. Atomic widths are due to the EiV — A7V 
conversion which in high angular momentum states is dominated by the one pion exchange. A joint 
analysis of the E~p — > An scattering data and atomic widths allows to extract the pion-hyperon 
coupling constant /j? AS /47r = 0.048 ± 0.005 (statistic) ±0.004 (systematic) or G^ AE = 13.3 ± 1.4 
(statistic) ±1.1 (systematic). 



I. INTRODUCTION 



It has been known for some time that hadronic atoms allow to test the long range part of hadron-nucleon interactions, 
the one-pion exchange (OPE), Ideal conditions are found in high angular momentum atomic states. There, 

the centrifugal barrier prevents the hadron to get close to the nucleus and the OPE effect is strongly enhanced over 
1^ the effects of two-pion or heavier meson exchanges. Detailed investigations are rather involved. In the best studied 
cases of antiprotonic atoms most of the long range pionic effects disappear due to spin and isospin averaging over 
the nuclear states of the bound nucleon. Decays of E-hyperonic atoms offer here a better chance, even though the 
\ experimental data are in general rather uncertain [p]-[5| . One exceptional case is the relatively recent measurement of 
Powers et al. Jj| that provides a precise width of 17(3) eV for the £~Pb atom in an "upper" n = 10, 1 = 9 state. It 
gives a good opportunity to study the role of one-pion exchange. 
The lifetimes of £ atoms arc determined by the rate of the reaction 

Y,- +p-> k + n (1) 

occurring at the nuclear surface. This process of hyperon conversion is dominated by the OPE and we exploit this 
feature to do a best estimate of the /^a_e coupling constant. Reaction (Q) has been studied in scattering experiments 
at low energies but the data are sparse and not very precise Differential cross sections are very uncertain. 

. £h ' These, as well as the elastic and other inelastic hyperon interactions, have been described by meson-exchange models 
|j3|--|l"Tl| . In addition to OPE, other short-range interactions are found to contribute. On the other hand, in high 
angular momentum atomic states of S atoms, the effects of short-range forces are strongly suppressed. The £ — N 
interactions in these atomic states are somehow equivalent to interactions in high partial waves of the free E— N system. 
However, there is a price for this simplicity: the atomic system involves the rather intricate nuclear structure. The 
latter determines the initial state of the atom and affects the final decay states. Fortunately, the conversion happens 
on a proton and the nuclear proton-density distributions are well known from \i atoms and electron scattering data. 
In addition, the proton momenta are quite well determined by nuclear shell models. The final states of the nucleus, A 
and neutron are less certain, however. These enter into the formalism, essentially in terms of the energy conservation 
that determines the energy and momentum transferred by the exchanged pion. In this paper it is shown that the 
atomic-level widths are very stable against these uncertainties. 

In section II of this paper we present the formalism and the determination of the 7rA£ coupling constant from 
the atomic widths. In addition to the Pb datum, the most accurate older results are also included into the best fit 
procedure. Some details of the OPE are given in Appendix A while Appendix B discusses the atomic wave functions. 
In section III limitations and uncertainties of our method are discussed and the error of the extracted 7rA£ coupling 
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constant is evaluated. These limitations stem from two sources: first the hyperon conversion process is not described 
by the pion exchange alone, there exist other amplitudes due to short ranged processes. Several models of these 
additional amplitudes are used and tested against the scattering data. The effect of such amplitudes on the atomic 
widths is found to be a 10% correction. The determination of the 7rAE coupling is weakly affected. Other uncertainties 
are due to the nuclear physics involved in the process. All these generate systematic errors which are comparable to 
the statistical one. 



II. THE NUCLEAR CONVERSION OF ATOMIC E HYPERONS 



The E~ hyperon, bound into atomic orbits, cascades down to states where it becomes converted into A. The 
quantum numbers and lifetimes of such states have been determined in X-ray experiments [|| , j| , || . ft turns out 
that these are levels of high angular momenta which force the conversion process to occur at the nuclear surface. The 
overlaps of the nuclear and atomic densities are localized in extreme surface regions and extend to distances as large as 
twice the nuclear radius Jl2| . Nuclear densities involved in the " upper" levels (states of the highest available angular 
momenta) amount to 3% of the central density. Processes of such peripherality allow for the standard low density 
simplifications: quasi-free scattering and single-particle picture of the nucleus. The final states may be described as 
free states. On the other hand, there exist a difficulty inherent to the surface studies related to the sensitivity to all 
ranges. For the problem of interest, we turn it into an advantage as it will allow us to precise the 7rAE coupling. 

This section presents the description of the E — A conversion mechanism. It follows the calculation performed in 
Ref. |l3[] for the antiprotonic atoms. First, a simple phenomenological picture based on the optical potential model is 
presented. Next, we elaborate a unitarity relation for reaction (|l|) and finally some technical details of this calculation 
are discussed. 

The tool to describe the atomic- level shifts and widths due to nuclear interactions is an optical potential, V opt . The 
simplest one is usually |12|] assumed to have the form 



yo P t {K) = f du p(R _ u ) ts7v(u) 

Meat J 



(2) 



where /zeat is the EAT reduced mass, p(R) is the nuclear density at a radius R and t-^jsr is the elastic scattering matrix. 
In Eq.(^|) the finite range in the T,N interactions induces a folding of the nuclear density p over this range. The atomic 
level widths are given by 

r/2 = - J dR ImV opt (R) | * E (R) | 2 

= — — [ dR [ du p(R-u) Im t SN (u) | tf s (R) | 2 (3) 
fJ-SN J J 

where \&e(R) is the atomic wave function that involves effects of nuclear interactions. Since ~ R and since 
the angular momenta I are large, the nuclear and atomic densities entering into Eq.(||) are well separated and it is 
only the longest range components of Im isiv( u ) that matter. To simplify the discussion we consider the widths 
of "upper" atomic levels, only. There are three advantages of this choice : i) the most precise data, ii) the largest 
I and iii) the almost unperturbed hydrogen- like atomic wave functions. The optical potential given by Eq.(|J) has 
some phenomenological success Q but for the purpose of this study we need a more microscopic approach that 
discloses the strength and range involved in Im t-^N- This quantity is related by unitarity to the inelastic reactions 
E~p — > An, E°n. In addition some specific atomic conditions must be met : the initial hyperon and nucleon are 
bound and then the E°n channel is blocked The unitarity condition is also affected by the presence of the nucleus. 
The latter takes away some recoil energy and may induce final state interactions of the outgoing hadrons. These 
interactions are neglected since the process is an extremely peripheral one, there are no Coulomb interactions and the 
final particles share about 70 MeV kinetic energy. On the other hand an assessment of the nucleon binding, of the 
nuclear recoil and of the rearrangement energies is of significance and will be studied below. 



A. The atomic-level widths and the unitarity condition 

The aim of this section is to calculate the rate of the nuclear E — A conversion. This is done in two steps: 
1) An amplitude for the Ep — An conversion is given in the leading order by OPE. This amplitude is introduced 
into the nuclear transition matrix element. Other "background" amplitudes, tb, will be discussed later. 
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2) The A emission probabilities are calculated and summed over states of A, n and over states of the final nucleus. 
For an isolated X p pair this procedure would produce the inelastic cross section and, via the unitarity condition, the 
absorptive amplitude Im tsN- Here, the summation over final states extends over free A, n states and over single 
proton-hole nuclear states. This summation generates Im t-£N and in the nuclear case it is folded over nuclear and 
atomic wave functions. These wave functions arise in the form of mixed single particle density matrices. 

Assume that the E hyperon in an n-th atomic state of quantum numbers [n, I, m] and a proton in a single particle 
state a convert into a free A and a neutron. The transition amplitude for this process is 

A n , a ,p, q = J dxdy * s (x)^(y) K(x-y) exp(iP.R + iq.r) (4) 

where ipf] (y) are the nuclear wave functions and V„ (x — y) is the OPE potential the expression of which is derived 
in Appendix A. The momenta and coordinates (P,R) refer to the center of mass while (q, r) denote the relative 
coordinates for the final pair. The final nucleus is specified, in the spirit of the impulse approximation, to be the 
initial nucleus left with a hole in the single particle state a. 

To calculate the absorption widths, amplitudes A n ^ a _p <q of Eq.Q) are squared, summed over the final nuclear states 
and integrated over the phase space of the A,n pair. In this way one obtains the following expression for the atomic 
level widths: 



T/2 = Y J J d(RR'rr') d(Pq)* s (R' + or') * S (R + or) <^,(R' - br') <p%(R - br) 



xK(r)V„(r') ^M-Eip-EM) exp( . R(R _ R/) + . q(r _^ (g) 

The coefficients a = Mn/{Mn + My) and b = Ms/ (Mjv + Ms) come from the transformation to the CM. coordinates. 
The energy excess Q = Ms — M\ — Eb{q) is due to the hyperon mass difference reduced by the nucleon binding 
Es{a). The values of Q fall in the range of 70 — 80 MeV, sizable by the nuclear standards and lead to simplifications 
that are valid at the nuclear surface. With typical £ and p momenta that are met in this region one finds the final 
CM. recoil energy E(P) to be about 10% of the relative kinetic energy E(q). Hence an average value for this recoil 
energy is used in Eq.(||). This allows to do the integration over P which generates a <5(R — R') and brings some 
simplicity to the otherwise involved integral. The calculations that follow are straightforward although the numerics 
has to be performed with care. 

The energy conservation fixes the final state relative momentum q which is also the momentum carried by the it 
meson. The actual value of q depends on the nucleon binding and recoil energies. Its central value is about 260 MeV 
but some 10% changes follow from the distribution of Eb{ol) + E(P) which characterizes proton states localized at 
the extreme nuclear surface. These changes are significant. 



B. Range effects 

In the limit of zero range EiV interactions, formula (^) for the conversion width may be expressed in terms of 
nuclear densities. For finite-range interactions the single-particle wave functions generate mixed densities. However, 
for simplicity and reference to the experimentally tested charge densities, one wants an expression in terms of the true 
densities. At the nuclear surface this can be done with a good precision, at least for the absorption rates summed 
over all nucleon states. The relation that allows it has been obtained in Ref. JTJ] and reads 

5^(R- bv') ip%(R- br) « p(R- bu) J nuc {k bw) (6) 

a 

where w = r r', u = (r + r')/2, p(R — bu) is the proton density. The J nuc (x) = Zj\{x)/x is an analog of Wigner 
function which describes the proton momentum distribution within the nucleus. The k is an effective local momentum. 
It is calculable in a shell model and at the nuclear surface it can be related to the proton separation energy [^5|. In 
a similar way we express the angular averaged atomic wave functions ^f^f by 

^-!— -^^(R + ar') * s (R + ar)«|* s (R + au) | 2 J at (aw), (7) 

m 

and J a t{aw) is obtained by an expansion of the atomic wave functions in terms of aw (see Appendix B). Collecting 
all formulas from (||) to (Q) one arrives at an expression 
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(v^/ttAw/ttAe) 2 UNA ( 1 _*4 2 „ 4 



x y rfRdu I * E (R + ou) I 2 p(R - bu) F fo i d (u). (8) 

The function Ffoid that describes the range of atomic-nuclear folding is given by 

, /" dflq exp(-m* | u + w/2 |) exp(-m* I u - w/2 I) 

Ffoid(u) = / aw- ; ■ 7 — j : j—. 

J 47r |u + w/2| I u — w/2 I 

x exp(iq.w) J nuc (bw) J at {aw). (9) 
To obtain the result of Eq. (@), few simple manipulations were carried out : 

1) As shown in Appendix A, the charged pion exchange involves some energy transfer that modifies the pion mass 
m w into to* . The Nf — 1.02 factor describes off-s hell effects for the bound nucleon. The pion exchange potential of 



Eq. (|A5|) , corrected for the S(r) singularity by Eq.(A6) has been used. The derivative operators are transformed by 
partial integration to derivatives over the initial and final state wave functions. These derivatives contain one large 
term due to the final An relative momentum q and smaller contributions due to relative momenta of £ and p. The 
large momentum produces a dominant g 4 term in the square bracket of Eq.(^). The effect of the initial momenta is 
obtained by calculating the derivatives over nuclear densities, atomic densities and correlation functions J nuc @ and 
J a t (0)- To the leading order in u/R, w/R these induce a change of q 2 into, 

q* 2 = q 2 + [a(l/R - l/Bn) + b/(2a p )f p } 2 ~ (kb/h) 2 , 

implemented in Eq.(|^). It amounts to a 25% correction, to obtain it we used the atomic wave function of Eq.( [B2| ) 
and a two parameter Fermi (2pF) nuclear density profile 1/(1 + e), where e = exp[(c — R)/a p ], c is the half density 
radius and a p is the surface thickness. In terms of the latter f p = e/(l + e). Other corrections due to J a t and terms 
(u/R) 2 , (w/R) 2 are found to be negligible. 

2) The spins of initial and intermediate baryons are averaged and summed. These generate the 1/3 and 2/3 weights 
in the square-bracket term of Eq.(||). 

Before performing the numerical calculations of Eq. (^)J) , let us elucidate the range involved in the folding function 
Ffoid- For the sake of argument let J nuc = J a t = 1 for these slowly varying functions. By going to Fourier transforms 
one obtains: 



16 f dfl q exp(i2u.p) 
rfoldW — — dp — — - — — — 

7r J Att (m* A + (p — q)- 2 ) (m u + (p- 



(10) 



where m! should be equal to m* but we have allowed for a more general situation which is met when several mesons 
of different masses are exchanged. From Eq.(|l0|) one infers that if m* = ml the integrand dependence on the (p.q) 



angle is rather moderate. The main bulk of Ff id(u) is given by 16n exp (— 1u\J to* 2 + q 2 J / \J m* 2 + q 2 . Its range is 

related both to the mass of the exchanged meson and to the momentum transfer, as expected from the uncertainty 
principle. More detailed calculations indicate some longer range contributions that induce small oscillations for very 
large u. For m' > to* formula ( |l0| ) generates an oscillatory behavior of Ffoid- This may be seen in the limit of large 
to' which yields 

„ sin(2qit) exp(— 2m*u) 

Ffoid = 167T ^P-^ ~^- 2 '-. (11) 

Zqu m z u 

Rapid oscillations in Ff D ia- at large values of u suppress such contributions to the atomic width. The m! ^S> to* terms 
are due to interferences of the pion exchange amplitude and other amplitudes of shorter ranges in space. These may 
arise from : exchange of heavier mesons, multiple scattering effects that are supposed to be of ranges shorter than 



1/2to w and vertex functions as shown in Eq.(A7) 



N ow, we attempt to describe the atomic level widths with the simple OPE mechanism including the form factor of 
Eq. dA7|). The most precise X-ray measurements of the upper level widths are collected in Table I. These results are 
reproduced by Eq.(||) with the best fit value / 2 ae/47t = 0.050(6) on the confidence level of xjdata = 2.38/4. In fact 
the fitted parameter is /irAfArZ-TrASj and to obtain the pion-hyperon coupling an average value of f 2 ^^/^ = 0.08 was 
used. These results were obtained with proton densities extracted from the 2pF charge densities of Fricke et.al |[L7) . 
The choice of these densities and the method to obtain the proton densities is discussed in more detail in the next 
section. The error of /^ax; given above is due to the experimental uncertainties. In addition, there exist uncertainties 
of the nuclear structure and simplifications in the conversion mechanism which induce systematic errors. To obtain 
these we analyze the T,p inelastic scattering data and discuss the basic nuclear parameters. 
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III. RESULTS 



The nuclear parameters of importance in this calculation are the energy excess Q and the nuclear recoil energy 
E(P). The energy excess is determined essentially by the separation energy of the valence protons. The recoil energy 
is given by the T,p pair momentum distribution, which follows the Fourier transforms of v E's(x) x (p^fx). These 
functions are localized around a point R a at the nuclear surface. Thus the Fourier transform is peaked at P w L/R a 
where L is the angular momentum of the pair with respect to the nucleus. The L is a sum of the atomic and nuclear 
single particle angular momenta L = l a + 1„. Since l a > l n the central value of L is l a and the most likely momentum 
is P = l a /R ~ l/m -1 . This estimates the recoil energy to be about 10 MeV. More detailed calculations indicate 
the spread of the recoil energies of some 5 MeV half-width. As discussed in the previous section an average energy 
D = (Eb{ol) + E(P)} is used to determine q, the momentum transferred by the pion, and the range involved in the 
folding function Ff id(u). The shell model generates rather small values for D and D/E(q) < 0.2. Nevertheless, 
uncertainties may arise, due to the way one averages over the binding energies Eb{o) or one calculates the recoil 
momenta P. These may be reflected in a few percent uncertainty in q. This problem is fortunately removed by a 
remarkable numerical stability of the atomic level width with respect to changes of D. This effect is related to the form 
of the pion coupling in Eq.(||). For an increased (decreased) q one finds a decrease (increase) in the atomic-nucleus 
overlap integral which is perfectly balanced by the increase (decrease) of q x (q* 4 ) in Eq.(||). This property holds on a 
few percent level for D in the region of 5 to 25 MeV. In this calculation we use an average value of D[MeV] = E s + 10 
where E s is the proton separation energy. Changing D by ± 5 MeV one obtains a change of 0.001 for the best fit 
/ttAE- This change is used as a partial estimate of our systematic error. 



A. Sensitivity study to short range amplitudes 

So far the Ep — An conversion has been described by the long range OPE amplitude. The justification of this 
approximation is the high angular momentum barrier that separates the atomic hyperon and the nuclear proton. An 
additional finding of the previous section is the suppression of an interference of this pion amplitude with amplitudes 
related to a shorter spatial extent (see Eq. ( |TT| ) . A check for these additional " background" amplitudes is possible since 
the free space conversion cross sections are known. In principle, one needs a full interaction model to do that. However 
the data on the YN interactions are sparse and the models || - |ll| differ in their content and their results. Here, 
we take a phenomenological attitude. A plausible form of the background amplitude, motivated by these theoretical 
models, is assumed. One or two free parameters are used to describe its magnitude. Next these strength parameters 
are fitted to reproduce the experimental conversion cross sections. The effect of this "background" amplitude on the 
atomic level width is calculated and the best fit to both atomic and scattering data is used to extract /^a^- Several 
forms of the background amplitude, tj, are used and the coupling constants obtained in this way are compared and 
"averaged". Such a procedure allows for an estimate of a systematic error in the coupling constant. 

The choice for additional amplitudes is based on several observations : 

(1) On the ground of Lippmann Schwinger equation one expects the scattering t(r) matrix to be of the form 

t(r) = V(r) + V(r) J dr' G(r - r') i(r') (12) 

where t(r) = t^(r) + th(r). The second term of this equation is due to multiple scattering effects. The "background" 
amplitude of interest may be due to the short range part of potential V(r) and/or to the multiple scattering contri- 
bution. What is significant is that the range involved in both these terms is expected to be at least shorter than that 
of the two-pion exchange amplitude. 

(2) The low momentum (140 — 170 MeV/c) inelastic E~p cross section of Engelman II shows no structure and no 
enhancement at low energies that would indicate a large scattering length in this system. The K- matrix parameteri- 
zation for a coupled AN, EA system in this state indicates the dominance of inelastic processes However, to be 
on a safe side, we fit the Massachusetts data (7l taken at higher momenta (200 — 580 MeV/c, 8 data points) to be 
more independent of the initial state interactions. 

We consider several models with different background amplitudes. The results of the best fits to the atomic and 
scattering data are summarized in Table II and plotted in Fig. |l|. The content of these models is as follows: 

A) OPE approximation is used to describe both the atomic widths and the reaction cross sections. These quantities 
are dominated by the tensor part of OPE. The atomic data by itself favor somewhat stronger /^s coupling, and this 
is reflected in a small reduction of the best fit in comparison to the result of the previous chapter . The overall fit is 
very good, and this leaves very little room for other contributions to the E — A conversion. 
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B) The two pion effects are simuf ated by Yukawa potential of the 1 / (2m„) range and f 2 coupling. It is supplemented 
with a form factor of source size parameter A = 1100 MeV. This force comes predominantly from a repetition of the 
tensorial force (see Eq. |l^) and contributes mostly to the spin triplet state, so the spin dependence of this potential 
is chosen to be (3 + a n -&y) / ^- The effect on the best fit f 2 ^ is minute. Similar results are obtained with an 
assumption of a tensorial form of this two pion force. The background amplitudes contribute some 10% to the atomic 
widths. A remark could be added at this point. Here one has a charged particle exchange and one knows in the NN 
sector that the uncorrelated charged two-pion exchange gives a small contribution, there it is dominated by the p. 
One could expect here a similar situation. 



C) A tensor component of the force (A4) is used with a large m! = m p mass of 770 MeV and a coupling f 2 



b ■ 



D) From SU(3) symmetry one exp ects a strong coupling with K meson as found in Refs. ||,fLl]l- The pseudoscalar 



exchange potential of the type (A5) is used with form factors A = 1100 MeV. Our best fit coupling constant f 2 = 
V^/knj: Jkna = -0.19 to be compared to the value - 0.387 of Ref. g and - 0.416 of Ref. @, although a much 
smaller number — 0.030 is obtained in Ref. Q. However, the "background" constants obtained here are just indicative. 
The data sample is not large enough for a real determination. These couplings contribute very little to the atomic 
widths and the corresponding \ 2 minima are quite shallow. The "background" amplitudes are used only to estimate 
the effects of possible short range forces compared to the basic longest range OPE. 

E) The p meson exchange amplitude of coupling 

(er x q).(<r xq) = i (S T q 2 + 1n N .a Y q 2 ) - (-S T q 2 - 2a N .a Y m 2 ) (13) 

is used with the mass m p = 770 MeV and the cutoff range A = 1100 MeV. 

F) The p + K meson exchange amplitude. This is a three parameter fit but only a marginal improvement has 
been obtained. Again the background amplitudes contribute about 10% to the dominant OPE. Now /j^ = —0.20 and 
ft I- 

These results show that the best fit value of f 2 ^ depends rather weakly on the uncertain short range interaction 
mechanism as can be seen in Fig. 0. There the errors shown for each model correspond to the quadratic sum of the 
statistical and systematic uncertainties discussed previously. From Table II the systematic error due to the imprecise 
knowledge of the short range interaction is found to be 0.003. 



B. The choice of nuclear densities 



Nuclear factors of prime importance are the proton density distributions p(r). These arc based on muonic atoms 
and electron scattering data. It is known that these experiments determine precisely R ms -the mean square radii 
of the charge distribution. However, for the problem of interest one needs higher moments of the proton density 
distribution, and these are not known very precisely. Since the atomic wave function behaves as r l one could expect 
the r 2 '-th moment to play the dominant role. More realistically, due to the exponential term in the wave function 



of Eq.(B2) it is the r 2 ' _2 -th moment of the density that gives the main contribution. According to Eq.(||) we are 
interested in the moments of proton distributions folded over a form-factor related to the pion exchange. Let us now 
discuss the relation of these quantities and the uncertainties involved. To be specific we study the best known nucleus 
40 Ca. 

The relation of bare (r 21 ) and folded (R 21 ) density moments may be obtained via studies of the corresponding 
Fourier transforms as given in Ref. p6[ . For I — 5 one obtains 

(R 10 ) = (r 10 ) + ^((r 8 )(u 2 ) + (r 2 )(u 8 )) +66«r 6 )(u 4 ) + (r 4 )(u 6 )) + (u 10 ) (14) 

where (u) are the moments of the folding function Ff id- Similar relations allow to express the proton density moments 



r 



2//\ 



in terms of protonic (u pc harge) and nuclear (R c harge) electromagnetic moments. The charge density moments are 
calculated from the experimental charge densities taken from the compilations of Refs. |17t| . One finds five different 
densities: two sets of 2pF and three sets of 3 parameter Fermi (3pF) profiles (these 3pF will be defined in the following 
paragraph). Next, we calculate the average proton density moments generated by those densities. These moments 
contribute to the folded moment of Eq. ( |l4| ) with the weights of 1, 0.71, 0.02, 0.33, 0.10 and 0.001 respectively. 
Thus the highest (r 10 ) moment contributes less than half of the whole folded moment. On the basis of these five 
experimental densities one can also calculate an error (r.m.s deviation) of each proton density moment. Thus the 
error of (r 8 ) is 11% while the error of (r 10 ) is as large as 29%. A total error of (i? 10 ) obtained from Eq. ( |l4| ) becomes 
18%. The same analysis done for the folded (i? 8 ) moment yields an error of 7.8%. These large errors offer certain 
upper limits since the errors of the (r 2fc ) are correlated and should not be added independently. 
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A different procedure to calculate the uncertainty is adopted here: for each experimental charge distribution we 
calculate a corresponding 3pF density. The parameters of this proton density are fixed to reproduce the lowest 
(H-charge)' (-^charge) aim (^-charge) moments. Next the atomic level width is calculated for each folded proton density. 
The mean square error of the average width obtained in this way amounts to 3.4%, only. With the same procedure 
we obtain errors of the other atomic level widths: 6.4% in Pb (two 2pF and two Gaussian charge density profiles), 
10% in Al (three 2pF profiles) and 11% in Si (two 2pF and two 3pF profiles). The overall uncertainty, weighted 
by the experimental errors, is estimated to be about 5%. For each nucleus we chose the proton density closest to 
the average and this density is used for further calculations. It turns out that in each case these are the proton 
densities based on the experimental charge densities by Fricke et.al (jl7|. The parameters (c, a p , w) found for the 
3pF ((1 + w(R/c) 2 /(l + e)) proton density profiles are: (3.01, 0.493, -0.01) for Al, (3.13, 0.479, 0.016) for Si, (3.69, 
0.518, -0.09) for Ca and (6.656, 0.474, -0.01) for Pb. 

The E hyperon in atomic orbit is affected by the complex optical potential of the nucleus. This changes the atomic 
wave function in the nuclear region. Two opposite effects contribute: i) the nuclear absorption described by Im V opt 
tends to suppress the wave function, ii) the Re V opt is weakly attractive at large distances and tends to enhance the 
wave function. Typical optical potential strengths are described by a scattering length bo- Phenomenological, best fit 
values are given in Ref. |12| and the simplest choice is bo = (.25 + i.l5)/m which characterizes an optical potential 
with the charge density profile. Solving Schrodinger equation for the level width in Pb one finds that the perturbative 
calculation underestimates the true width by about 11%. The effect of attraction prevails over the repulsive effect 
due to absorption. The same exists also within our description of the absorption. However, the effect is smaller since 
our absorptive pion exchange potential is longer ranged that the real potentials. 

To account for the multiple scattering effect we express Eq.(||) as the perturbative result for a folded potential as 
indicated by Eq.(|^). A shift of an argument by cm is necessary to perform this. Next, such a potential is supplemented 
by the real optical potential discussed above. The level widths obtained from the Shrodingcr equation are larger than 
the perturbative level widths by 5.8% in Pb , 2.9% in Ca , 3.3% in Si and 3.1% in Al. The net result of this 
calculation is that the coupling constant / 2 AS extracted from the widths becomes smaller than the constant extracted 
in a perturbative way. The difference amounts to 4%. These corrections have been included in the final results given 
in Table II. 

The optical potential correction procedure is not unique as a number of real optical potentials exist in the litera- 
ture [ p2| . Different values of bo induce an additional 1% uncertainty. This we lump together with the uncertainty due 
to the charge distribution . Now , the total amounts to a 6% effect. 

Thus, nuclear densities induce the main part of our systematic error. The corresponding error of 0.003 arises in the 
•/ttAE coupling constant. 

The mixed nucleon density matrix of Eq. (^) introduces no serious uncertainties. As the ranges involved are given 
by 1/ ^2-\/q* 2 + to* 2 ^ one needs this expression for w < 1.5 fm. The whole effect of J nuc amounts to a reduction of 
r by less than 5% and the uncertainties in Eq.(K3|) are negligible. 



IV. CONCLUSIONS 



It has been shown that the upper level widths in E hyperonic atoms present a good source of information on the 
/wAE coupling constants. The complications due to nuclear physics may be kept under a reasonable control. The 
experimental uncertainties still present a problem. 

Averaging the best fits gives a value of 

/^As/ 47r = °- 048 ± °- 005 (statistic) ± 0.004 (systematic) . (15) 

This corresponds to G 2 AS /47r = /| as [(Ma + Ms)/m w ] 2 /4tt, i.e. 

G^ as /4tt = 13.3 ± 1.4 (statistic) ± 1.1 (systematic) . (16) 

This result is obtained here from the combined analysis of the atomic widths and the inelastic hyperon scattering cross 
sections. This result differs from the numbers used, or obtained in the analyzes of the scattering data. In particular : 
/J AE /47T = 0.041 has been used by Nijmegen group (|J^] while the Julich group Q obtains / 2 as /4tt = 0.034 and/or 
/ 2 AS /47r = 0.028 It turns out that the atomic data favor larger values of the pion hyperon coupling. Had we used 
a lower value for the pion-nucleon coupling constant such as that determined by the Nijmegen group Q viz., 0.075, 
instead of the classical value 0.080, our /^ As /47r would even have been 7% larger. 

The errors of this determination, summarized in Table III, consist of four components: i) the statistical error due to 
the data uncertainties, 0.005, the systematic error due to ii) uncertainties in the proton density distributions, 0.003, 
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iii) uncertainties of the energy release in the capture process 0.001, iv) the uncertainty due to poor knowledge of 
the short range hyperon-nucleon interactions estimated to be about 0.003. Altogether, an error of 0.007 is obtained. 
In order to reduce the large statistical part, both new measurements of the £ atomic widths and free "E~p — ► An 
scattering data would be welcome. 



We thank D. Vauthcrin for helpful discussions and S. W. is indebted to J. Dabrowski for useful comments. The 
authors acknowledge support of IN2P3 /INS collaboration agreement and KBN Grant 2P03B 048 12. 



APPENDIX A: PION EXCHANGE FORCES 

To make the paper self-contained several useful formulas are collected in this appendix. These concern the pion 
exchange in the EAT — > AA~ conversion. Due to the Ms — Ma mass difference this process involves momentum and 
also some energy transfer. It makes the exchange force to be longer ranged as compared to the N — N interactions. 
The E~ and p being bound into the nucleus arc not on mass-shell Dirac spinors. To derive the OPE potential we shall 
nevertheless consider them as free on mass-shell spinors. The reduction of the pseudo-vector 757 M <9 M 7r pion-nucleon 
(ttN) and pion-hyperon (ttY) couplings leads, in the momentum representation, to 

t t I \ VZfirNN T/ JttAS t , 1 / a 1 \ 

VM) = Vn — Vy w^w^eytk (A1) 



where 



and 



V N = J (E n + M N )(E P + M N ) a N . ( 



M N E P + M N 



(A2) 



V Y = Mz + M A (E A + M N )(Ey + M N p__ 

2 V M a M s \Ea + M a E-s+Mx' 



Here p and p' denotes the initial and final CM. momenta of the reaction (111), E n = Wp' 2 + M N , E p = ^/p 2 + M'^, 



E A = J p' 2 + M 2 and E% — ^/p 2 + . In the free space AE = E n — E p = E% — E A . If one neglects p' , p 2 terms 



and A — S mass difference, Eq. (JAlj) with Eqs. (A2) and (A3) reduce the corresponding momentum representation 
OPE to 



m~ — (AE)' S ™ l 



where q = p' — p is the momentum transfer. In the atomic state AE = Ey,n — £ajv is the energy transferred by 
mesons. The AE of w 40 MeV at low energies makes an effective pion mass m* = \Jm^ ~ (AE) 2 of about 135 MeV. 
The coupling constant appropriate for S^p — > Ap reaction is f 2 = \/2 f^NN fir as- In principle, it is the hyperonic 
constant /„ae which is to be determined here. In reality, this determination is affected by a small uncertainty of the 
pion nucleon constant /^jv n Pq] ■ 



In the coordinate representation the potential of Eq. ( A4 ) corresponds to 



mn a 4z7T7* 



This potential is strongly singular at the origin, and this singularity has to be removed by a zero in the wave function 
induced by a r epulsive core potential. A simpler and more customary way is t o remove the S(r) like term from 



expression ( |A5| ), |19[. This is conveniently done in momentum representation (A4) by the splitting 



(<r A r.q)((Ty.q) = i (S T q 2 + o- N .er Y q 2 ) ^ (St q 2 ~ ct n .(Ty m* 2 ) (A6) 



8 



where St — 3(trj y.q) (<Ty.q) — a^.ay is the tensor force that does not induce the 6(r) singularity. The step indicated 
by arrow in Eq.(A6) describes the subtraction of 8(r) CT^.rry f /m„ from the spin-spin potential. For the S — A 
conversion one finds q > 2m* which makes the tensor force to be the strongly dominant one. 

The size of pion source is described by a vertex function F^b b = (A 2 — m* 2 )/ (A 2 + q 2 ) that induces a change in the 
pion potential V 7T {q) =>■ V K -(q)FT V BBF%A'S- The same source size parameter A = 770 MeV i s us ed for the nucleon and 
hyperon vertices. These vertex functions induce a change in the spatial form of potential (|Aq) 



exp(— m*r) exp(— m*r) 



A 2 - 



2A 



- I cxp(— Ar) 
r ' 



(A7) 



This finite source potential is weaker than the point source one. However, in atomic states of high I the effect of 
source size is very small as the centrifugal barrier suppresses short distances. 

Keeping only first order corrections in the atomic states one has E n ~ Mjv + p' 2 / (2Mn), E p ~ Mat — Eb, 
Ea ~ A/a + p' 2 /(2Af A ) and £ E ~ Af s . Eqs. (|aJ) and reduce to 



V N 



E B 
4M 



N 



8M N 



o-jv-q 



and 



V Y 



1 



M s - M a 
2(M S + M A ) 8Af A 



cry.q 



It can then be seen that the coupling / of Eq.(A5) is renormalized by 

/ Mv- Ma R d n' 2 

N f =[l 



2(M E + M A ) 4M N 8M A 8Mat 



(A8) 



£Jb and p' take slight different values depending on the nucleus considered but on the average iV/=1.02. The /jrAS 
values quoted in this work have been corrected by this renormalization. 

The cross section for > A + n within the Born one pion exchange approximation becomes 



da q\N ( NN fir 



dfl gsTV \ 4-7T m 2 

where q\N (<?EAf ) are the final(initial) cm. momenta. 



r irBB p ttAS 



m*73 + 2 g 4 /3 
(m* 2 + g 2 ) 2 



(A9) 



APPENDIX B: ATOMIC DENSITY MATRIX 



The unitarity expression for level widths requires some simple representation of the mixed atomic density 

p at (Z, z) = ^ J2 *(Z + »)*(Z - 



(Bl) 



In the text we have Z = R+a u and z = a w/2. The p a t(Z, z) is evaluated for small z/Z and the expansion parameter 
becomes e = z 2 /Z 2 w (b/2) 2 w 2 /R 2 . For Pb atom this quantity is very small, e < 0.002, but the atomic wave function 
change rapidly in space and some care is necessary. 
The states of interest are the circular orbits and 



,(Z) = N nl Y lm {Z) Z l exp( -) 

nB 



(B2) 



where B is the Bohr radius and N n i is a normalization. To perform the expansion let us notice that (Z + z).(Z — z) = 
Z 2 (1-e), |Z + z||Z-z| = Z 2 (1 + 26-46C 2 ) 1 / 2 and 

(Z + z).(Z-z) 2 

t= rr^z r - = X ~ 1 — 2e + 2ec 

Z+z Z-z 
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where c = cos(Zz). The most important correction in Eq. (Bl) is due to spherical harmonics as the hyperon momenta 
are predominantly tangential to the nuclear surface. One has 

A ~T- 

^-^Y, Y im(Z + z)Y lm {Z -z) = Ptix) « P t (l - 2e + 2ec 2 ) wl-l(l+ l)e(l - c 2 ) (B3) 

where we use P/(l) = 1(1 + l)/2. The other pieces of atomic wave functions expanded in a similar way produce much 
smaller corrections. Thus \Z + z\ l \Z - z\< w Z 21 (1 + Ze(l - 2c 2 )) and 

/ |Z + z| + |Z-z|\ / -2Z\ f e 
exp LJ « exp 1 



nB J \ nB J V nB 

Altogether one obtains 



p at (Z,x) »p«t(Z,0) Pi(l-2e + 2ec 2 ) (1 + le(l - 2c 2 )) 1 + — 



« p ot (Z, 0) - 4e/3) (1 + /e/3) (1 + —J (B4) 

where in the last line an average value of c 2 = 1/3 was used. In practice it is only the Pi(l — 4e/3) term that is the 
important one giving a few percent correction to the atomic density p at (Z,0). 
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TABLE I. Experimental and calculated results for upper level widths, r(eV), in circular orbits. I is the level angular 
momentum. 



Atom 


E 


k 


f. 1 


-Tcxperiment (b"V) 


-^calculated 


(eV) 


"Al 
28 Si 
40 Ca 

208p b 


4 
\ 

5 


4 
4 
5 
9 


0.24(6) 
0.41(10) 
0.41(22) 
17(3) 


0.21 
0.51 
0.17 
15 






TABLE II 


Summary of the best fits to the atomic and scattering 


data. 




Model 


A(tt) 




B(tt + Vi(2tt)) 


C(n + V T {p)) D(n + K) 


E(vr + p) 


F(tt + K + p) 


atomic %*/ data 
full x 2 /data 
fb 

/ttAE 


2.42/4 
5.90/12 


0.047(4) 




2.42/4 
5.84/12 

-0.09 
0.045(4) 


2.33/4 2.31/4 
5.74/12 5.67/12 

2.4 -0.19 
0.050(4) 0.051(5) 


2.43/4 
5.89/12 

0.4 
0.047(5) 


2.28/4 
5.58/12 
-0.20 (K); 1.7 (p) 
0.051(5) 



TABLE III. Summary of errors in the determination of /^ae- 



Source 


Energy 


Proton density 


Background 


Statistics 


of error 


release 


distribution 






Error 


0.001 


0.003 


0.003 


0.005 
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0.058 



0.053 



5+H 



« 0.048 



0.043 



0.038 
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FIG. 1. Results of the best fits to atomic and scattering data with different background amplitudes (see text). Errors are 
the quadratic sum of statistical and systematic uncertainties 
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